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Chapter  1 
INTRODUCTION 

1.  Forex^ard:   It  is  often  maintained  that  thermodynamics  Is 
concerned  with  reversible  processes  and  equilibrium  states,  and 
that  problems  associated  with  irreversibility  lie  outside  its 
scope.  Starting  with  three  fundamental  postulates  and  a  number 
of  state  functions  defined  for  the  equilibrium  states,  the 
classical  equilibrium  thermodynamics  has  developed  into  a  most 
satisfactory  axiomatic  science.  Another  aspect  of  the  equili- 
brium thermodynamics  is  its  phenomenological  attitude.  With  its 
foundation  seated  on  a  basis  free  from  any  considerations  of  the 
structure  of  matter,  the  equilibrium  thermodynamics  has  yielded 
important  results  entirely  unshaken  by  the  quantum  revolution. 
There  arises  the  question  of  the  feasibility  of  the  extension  of 
the  methods  of  equilibrium  thermodynamics  into  the  realm  of  non- 
equilibrium  (or  irreversible)  phenomena.  A  phenomenological 
theory  of  Irreversible  processes  should  aspire  to  encompass 
within  a  simple  scheme  a  wide  range  of  non-equilibrium  phenomena. 
Like  other  macroscopic  disciplines,  such  as  fluid  mechanics  and 
electromagnet Ism,  the  thermodynamics  of  Irreversible  processes 
should  be  based  on  a  scheme  of  continuum  field  with  characteris- 
tic continuous  functions  of  space  coordinates  and  time.   In  con- 
trast to  the  classical  equilibrium  thermodynamics,  In  which  the 
state  variables  are  usually  independent  of  the  space  coordinates, 
the  thermodynamics  of  Irreversible  processes  should  have  Its 
formulation  In  the  forms  of  local  equations  just  as  in  fluid 
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dynamics  and  the  Maxwell! an  theory  of  electromagnetisra. 

2.   On  Irreversibility :   There  are  two  erroneous  criteria 
popularly  considered  to  be  characteristic  of  Irreversibility, 
namely  spontaneous  occurrence  of  a  process,  and  the  running  of 
the  process  at  a  finite  rate.   Bridgman  (1)  has  Indicated  that 
a  finite  and  reversible  flow  of  heat  could  be  maintained  by  the 
action  of  a  thermo-motlve  force  (in  analogy  with  the  electro- 
motive force  of  electrodynamics).  The  increase  of  entropy  is 
the  ultimate  criterion  of  irreversibility.  There  are  probably 
other  criteria  of  less  general  scope.  An  explicit  criteria  for 
irreversibility  suggests  itself  when  the  differential  equations 
describing  some  time-dependent  physical  processes  are  examined. 
If  those  equations  are  invariant  with  regard  to  the  algebraic 
•  sign  affixed  to  the  variable  time,  the  process  is  considered  as 
reversible,  otherwise  it  is  called  irreversible.  As  an  example, 
the  wave  equation  which  describes  the  propagation  of  waves  in  a 
non-absorbing  medium  is  (2) : 

i   a2u  =  d2u    a2u  .  a2u  M  ,, 

7^7^    r?    7~z    7"?  (1"1) 

Cot    dx    3y    dZ 

It  is  seen  that  this  equation  is  invariant  for  the  substi- 
tution of  t  by  -t  and  the  propagation,  which  is  of  natural 
occurence  and  takes-  place  at  a  finite  rate,  is  a  reversible 
process. 

On  the  other  hand,  the  transport  equation  which  character- 
izes the  usual  physical  processes  of  diffusion  and  conduction  of 
heat,  etc., 


i  3t  B  a2T  ,  a2T  .  a2T  ,.  Q, 

-y -ir  +   — r  +  — y  (i-2) 

is  not  Invariant  with  respect  to  such  a  substitution  and  it 
describes  an  irreversible  process,  or  more  exactly,  the 
irreversible  approach  to  equilibrium. 

3.   Historical  Development:   There  began  to  appear  in  the 
middle  of  the  nineteenth  century  scattered  treatments  on  irre- 
versible processes.  The  study  on  the  action  of  thermocouple  by 
W,  Thomson,  Lord  Kelvin  (3)  was  probably  the  first  case  known  in 
the  scientific  literature.  Lord  Kelvin  suceeded  in  correctly 
deriving  two  relations  by  resorting  to  an  altogether  intuitive 
and  unjustifiable  inference,  known  as  Kelvin's  Hypothesis.   In 
the  action  of  thermocouple  there  arc  two  irreversible  processes 
accompanying  the  thermoelectric  process,  namely  thermal  conduc- 
tion and  the  development  of  Joulean  heat.  Kelvin  postulated  that 
reversible  process??  ueve   separable  from  irreversible  processes 
in  the  sense  that  thermodynamics  could  be  applied  to  reversible 
processes  as  if  the  accompanying  irreversible  processes  were 
absent.  In  Kelvin's  work,  the  thermoelectric  phenomena  was  taken 
as  the  reversible  process  to  which  application  of  thermodynamics 
was  directed.   It  may  be  cited  that  Boltzmann  (k)   attempted  with- 
out success  to  justify  Kelvin's  Hypothesis. 

It  was  not  until  193 1  when  Onsager  (5)  proposed  the 
"reciprocal  relations"  that  the  treatments  of  irreversible 
phenomena  were  founded  on  a  rational  scheme.   Later  on  Meixner 
(6)  and  somewhat  later  Prigogine  (7)  succeeded  in  establishing  a 


k 

consistent  phenomenological  theory  of  the  irreversible  processes, 
incorporating  both  Onsagerts  reciprocal  relations  and  a  scheme 
for  explicit  calculation  of  the  so-called  entropy  source  strength 
in  a  certain  number  of  physical  situations.   In  this  direction  a 
unified  front  was  rapidly  developed  in  the  last  two  decades. 

k.     Onsagerts  Reciprocal  Relations:  In  its  original  form, 
the  context  of  Onsagerts  theory  (£)  may  be  stated  as  follows: 

When  a  set  of  n  thermodynamic  fluxes  J.  are  linearly  related 
to  n  thermodynamical  forces  X, : 


n 


Ji   £  Lij  xj  ,     i  =  1,  2,  .  .  .,  n,  (1-3) 


and  if  each  flux  is  the  time  derivative  of  a  thermodynamic 
variable  a-_: 

Ji  ■   *L  (l4) 

while  each  a^  is  an  even  function  of  the  particle  velocities 
(for  instance,  local  energies,  concentrations,  etc.  which  are 
not  affected  by  the  sense  or  direction  of  velocity  vectors), 
and  the  entropy  of  the  system  could  be  approximated  by: 

,   n   n 
S  =  *  Z   Z  g.j  a.  a  (i„5) 

d   i=l  j=l   J   L  J 

and  the  forces  x.  have  been  chosen  so  that: 
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x.  =(2£4=  2  g.  .  a.  (1-6) 

then,  the  matrix  L. .   of  the  phenomenologlcal  coefficients 

L.  .  in  Equation  (1-3)  Is  symmetric;  i.e., 

L.  .  =  L..  (1-7) 

The  modification  required  In  the  general  case  (when  the  a. 
can  be  odd  functions  of  the  particle  velocities,  e.g.,  momentum 
densities)  were  given  by  Caslmler  (8). 

Under  certain  circumstances  the  matrix  L,j   Is  not  symme- 
trie.  When  an  external  magnetic  field  acts  on  the  system,  the 
reciprocal  relations  become: 

LLj  (3)  =  LjL  (-B)  (1-8) 

where  B  Is  the  magnetic  Induction.  Thus,  L. .  Is  the  same  func- 
tlon  of  B  as  L_. .  Is  of  -B. 

j  !•  » 

A  similar  relation  holds  for  the  situation  In  which  the 
Cor lolls  forces  are  present,  I.e.,  when  the  system  Is  In  rota- 
tion. The  reciprocal  relations  In  this  case  are: 

Lu  (do)   =  Ljl  ("w)  ("l-9i 

where  co   Is  the  angular  velocity  of  rotation. 

Some  elaboration  on  Equation  (1-5)  will  probably  clarify 
the  approximations  Involved.  when  the  entropy,  S,  of  a  system 
Is  considered  as  a  function  of  thermodynamic  variables  a.: 
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S  =  S(a1,   a2,    .    .    .,   an)  (1-10) 

where  all  the  quantities   represent  the   deviations  from  their 
equilibrium  values,   S°  and   a0    ,    i   =   1,    2,    .    .    .,    n. 

Taylor  series    expansion   of  S  around  a0,    leads    to 

ft  n       n  2 

1=1    W   a1=0        L       2  1=1  J-X   '^J3SJ)   ai=0,a.=0   (1-11) 

•  a,  a.  +    .    .    . 

where  terms  of  degree  higher  than  the  second  are  neglected. 
Since  the  entropy  of  the  system  assumes  the  maximum  value  at 
equilibrium  condition,  the  first  term  of  the  right  had  side  of 
Equation  (1-11)  vanishes,  and  Equation  (1-11)  is  reduced  to 
Equation  (1-5) 


1  n   n 

2  i=l  j=i  ^      L   J  (1_5) 


where 


V   :     a  =0,a.=0 


An  important  theory  connected  with  the  application  of 
Onsagens  method  will  be  stated.  Equation  (1-3)  represents  a 
mathematical  statement  of  the  physical  phenomena  to  be  called 
"coupling"  hereafter.  Two  rate  processes  are  said  to  be  coupled 
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when  the  force  of  one  of  them  Induces  flux  of  another.  Thus  a 
temperature  gradient  Ls  known  to  cause  diffusion  of  matter  (Soret 
effect)  and  on  the  other  hand  a  concentration  gradient  tend3  to 
induce  heat  flow  (Dufour  effect).   Equation  (1-3)  describes  the 
flux  of  a  certain  kind  (heat  or  material  of  certain  components) 
In  terms  of  linear  combinations  of  all  the  "forces"  present  in 
the  system.  The  coefficients  L. .  characterize  the  coupling 
effect.  The  vanishing  of  a  coefficient  L, .  indicates  the 
absence  of  coupling  effect  between  two  phenomena.  There  is  a 
rule  known  as  the  Curie's  theorem  (9)  which  offers  a  definite 
answer  as  to  the  occurrence  of  coupling  between  two  phenomena. 
Curie's  theorem  states  that  entities  whose  tensorial  characters 
differ  by  an  odd  integer  cannot  interact  in  an  isotrophic  media. 
Some  physical  processes  are  often  classified  according  to  their 
tensorial  character.  The  tensorial  character  of  a  process  is 
datarsalned  by  tensorial  rank  of  its  force.  Under  the  category 
of  the  scalar  phenomena  (tensor  of  rank  one)  there  are  chemical 
reactions  and  structure  relaxation.  The  vectorial  processes 
(tensor  of  rank  two)  include  diffusion,  heat  conduction,  etc. 
Viscous  phenomena  is  an  example  of  the  tensorial  (tensor  of  rank 
three)  processes. 

In  practice,  however,  the  coupling  effect  between  phenomena 
of  different  tensorial  ranks  is  always  negligibly  small.  Con- 
sequently, only  coupling  among  processes  with  equal  tensorial 
ranks  needs  to  be  considered. 

5.  Significance  of  Thermodynamic  State  Functions  in  Non- 
equilibrium  States  :  The  thermodynamics  of  irreversible  processes 
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aims  at  extending  the  method  and  scope  of  the  classical 
equilibrium  thermodynamics.  All  the  thermodynamical  state  func- 
tions (such  as  temperature,  entropy,  etc.)  which  are  strictly 
defined  for  equilibrium  states  are  to  be  used  in  the  description 
of  non-equilibrium  phenomena.  This  rather  contradictory  atti- 
tude may  be  justified  when  the  treatments  are  confined  to 
phenomena  which  are  not  far  removed  from  equilibrium  states. 
In  this  conjunction,  a  postulate  which  may  be  named  the  postu- 
late of  local  equilibrium  (9)  will  have  to  be  introduced.  This 
may  be  stated  as  follows:  for  a  system  in  which  irreversible 
processes  are  taking  place,  all  thermodynamic  functions  of  state 
exist  for  each  element  of  the  system.  These  thermodynamic 
quantities  for  the  non-equilibrium  system  are  the  same  functions 
of  the  local  state  variables  as  the  corresponding  equilibrium 
thermodynamic  quantities.  The  validity  of  this  postulate  has 
been  investigated  by  Prigogine  (10)  from  a  microscopic  viewpoint. 
In  the  Chapman-Ens kog  (11)  theory,  the  molecular  distribution 
function  is  expressed  in  an  infinite  series 

f  =  f0  +  fl  +  f2  +  '  '  •  '  (1-12) 

where  fQ   is  the  equilibrium  molecular  distribution  function,  f 
is  the  first-order  corrections,  and  f£,  etc.,  are  subsequent 
corrections.  For  a  system  which  is  close  to  the  equilibrium, 
tiie  successive  corrections  f^  f£,  etc.  become  smaller  and 
smaller.  Prigogine  (10)  has  established  that  the  usual  thermo- 
dynamic quantities  defined  for  the  equilibrium  conditions  retain 
their  significance  when  the  molecular  distribution  functions  can 
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be  characterized  by  f  =  f  +  f    Prigogine's  analysis  was  based 
on  mixtures  of  monoatoraic  gases  and  dependent  on  the  Chapman- 
Enskog  model.   In  general,  the  domain  of  validity  of  the  postu- 
late of  local  equilibrium  is  still  not  established. 
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Chapter  2 
ON  THE  APPLICATION  OP   ONSAGER»S   RECIPROCAL  RELATIONS 

There   is   a  widespread  misunderstanding  concerning  the   choice 
of  fluxes   and  forces    in  conjunction  wtth  the  use    of  the  Onsager 
reciprocal  relationships.     For  the  sake   of  convenience,    the  con- 
text  of  Onsager»s   relations,    as   stated   in  Chapter  1,  will  be 
reproduced   here. 

When  a  set   of  n  "thermodynamical  fluxes"   J.    are   linearly 
related  to  "n  thermodynamical  forces"  X, : 

n 
Ji  =  ^  Llj  XV  i  =  1,    2,    .    .    .,    n  (2-1) 

and  if  each  flux  is  the  time  derivative  of  a  thermodynamic 
variable  a.  : 

Ji  "  V  (2-2) 

while  each  aL  is  an  even  function  of  the  particle  velocities  and 
the  entropy  of  the  system  could  be  approximated  by: 

.     n       n 

s  =?i?i  i?i  su  Vr  (2-3) 


and   if  the  forces  X.   have  been  chosen  so  that:. 

J 

J         aai  i=l     LJ      *•  K     ^' 


• 

a 

s 

-  z 

v    T 

1          1- 

1=1 

11 

then,  the  Matrix  L. .   of  the  phenomenological  coefficients  L. . 
in  Equation  (2-1)  is  symmetric;  i.e., 

L.  .  -  L„  (2-5) 

A  celebrated  equation  relating  the  entropy  production  to  the 
"forces"  and  "fluxes", 

(2-6) 

may  be   obtained  as   follows.      One   differentiates  Equation   (2-3) 
with  respect   to  time   to  obtain 

n       n 

S  =    Z     Z    g,,  a  a.  (2-7) 

1=1   j=i      1J      I    3 

On  comparing  this  expression  with  Equations  (2-2)  and  (2-lj.),  one 
sees  that  Equation  (2-6)  is  a  necessary  consequence. 

Although  Equation  (2-6)  results  from  Equations  (2-2)  through 
(2-Ii_),  there  is  no  assurance  that  the  Onsager  Reciprocal  Rela- 
tions as  expressed  in  Equation  (2-5)  follow  from  Equations  (2-1) 
and  (2-6).   In  current  literature  (12),  (13),  one  encounters 
many  examples  of  indiscriminate  choice  of  forces  and  fluxes 
based  on  the  erroneous  notion  that  Equations  (2-1)  and  (2-6) 
establish  Onsager* s  Reciprocal  Relations.  Coleman  and  Truesdell 
(11+)  have  demonstrated  that  a  set  of  conjugated  forces  and  fluxes 
satisfying  Equations  (2-1)  and  (2-6)  may  not  fulfill  the  require- 
ment of  Onsager* s  Reciprocal  Relations.   In  the  following,  a 
brief  proof  essentially  after  Coleman  and  Truesdell  will  be 
given.  Based  on  a  set  of  conjugated  fluxes,  J,  and  forces, 
Xj,  which  obey  Equations  (2-1),  (2-6),  and  (2-5),  one  can 
construct  in  an  Infinite  number  of  ways  new  sets  of  fluxes  and 
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forces  satisfying  Equations  (2-1)  and  (2-6) : 


J,  -  £  L./  X  (2-10 

1    i=l  LJ    0 


n    f   ,    n 

Z  X.   J.   =  Z  X.J.  =  S  (2-6«) 

1=1  l   L    1=1  l   L 

but  yielding  nonsymmetric  matrices  [l. .   of  the  phenomenological 


coefficients  L  . ..  To  be  specific,  let  f W.  ,]  be  any  nonzero 
n  x  n  skew  matrix;  i.e.,  let  the  components  of  hi.  .]  be  such 
that  for  all  the  pairs  of  i,  j: 

V   s  -W   •   1  =  1   P  w 

(2-8) 
j  =  1,  2,  '.  .  .,  H 

and,    for  at   least   one  such  pair  Ttfith   j  ^   1,    let  W,  .  ^   o. 

Consider  the   set   of  fluxes   J.    and  forces  X.    defined  as 
follows : 


<sV|iVxJ  (2'9) 

Xi  =  XL  (2-10) 

for  all  i,  j.   Then: 


n   |   |    n         n   n 

Z i  J.  x,  ■  Z    J.x  +  Z    Z  v. ,  x,  x.  (2-u) 

1=1  L  l       1=1  *•  L   1=1  1=1  5- J   J   L 
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That  the  second  terra  of  the   right  side   of  the  above  equation 
vanishes    Is  a  consequence   of  Equation   (2-8).      Hence,    it  follows 


that: 


a       ,      ,  a 


Z     JL  X,    =    Z     JL  X-    =  S  (2-12) 

1=1      l      X        i=l     L      X 

Ifow,  from  Equations  (2-1),  (2-9),  and  (2-10),  there  results: 

and  one  obtains,  by  comparing  Equation  (2-13)  with  Equation 
(2-11) 

Lij  =  Lij  +  ¥ij  (2-Hi) 


or 


IV I  "  fLij)  +  [WUJ  (2-15) 


t 


Therefore,   L.  .   ,  being  the  sura  of  a  symmetric  matrix  L,  . 
«  ij 

and  a  nonzero  skew  matrix  w. .  ,  is  nonsymmetric. 

Host  of  the  current  treatments  (1|) ,  (12),  (13),  on  the 
application  of  Onsager's  theory  rely  on  Equation  (2-6)  in  the 
choice  of  conjugated  forces  and  fluxes. 

It  has  been  recognized  that  phenomenological  theories  of 
various  kinds  easily  yield  a  relation  of  the  type  Equation  (2-6) 
while  relations  of  the  type  exemplified  in  Equations  (2-2) 
through  (2-4)  are  more  difficult  to  establish.  Therefore,  if  one 


Ill- 
takes  Equations  (2-1)  through  (2-1'.)  as  the  hypothesis  of 
Onsager's  theorem,  then  the  choice  of  the  fluxes  and  forces  is 
greatly  limited.  Of  course,  it  is  possible  that  some  other 
choice  of  fluxes  might  also  lead  to  Equation  (2-5).  But,  until 
the  forces  and  fluxes  in  any  particular  case  are  physically 
identified  "by  some  property  more  specific  than  their  mere  entry 
into  the  bilinear  form  of  Equation  (2-6),  application  of 
Onsager»s  theorem  is  precarious. 
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Chapter  3 
ISOTHERMAL  MULTICOMPONENT  DIFFUSION 

1.   Introduction:  The  main  context  of  the  previous  chapter 
was  to  demonstrate  that  the  legitimate  choice  of  conjugated 
forces  and  fluxes  in  connection  with  Onsager's  Reciprocal  Rela- 
tions rests  upon  a  more  fundamental  criterion  than  the  invari- 
ance  of  entropy  production.  Use  of  the  set  of  original 
conditions  as  expressed  by  Equations  (2-2)  through  (2-Ij.)  entails 
a  detailed  knowledge  of  microscopic  nature  and  will  defy  the 
spirit  of  phenomenological  approach.  Prom  a  phenomenological 
point  of  view,  the  correct  choice  of  conjugated  forces  and 
fluxes  relies  solely  upon  experimental  investigations.  A  compre- 
hensive experimental  research  in  an  attempt  to  examine  the 
validity  of  Onsager's  Reciprocal  Relations  has  been  carried  out 
recently  by  D.  G.  Miller  (13).   In  his  experimental  investigation 
Miller  selected  the  conjugated  forces  and  fluxes  on  the  basis  of 
the  invariance  of  entropy  production.  As  was  pointed  out  in 
Chapter  2,  the  criterion  of  the  invariance  of  entropy  production 
as  expressed  in  Equation  (2-6)  is  an  erroneous  one.   In  this 
light  one  may  look  at  Miller's  work  as  a  search  for  legitimate 
choices  of  conjugated  forces  and  fluxes  rather  than  a  direct 
attempt  at  the  verification  of  the  Onsager's  Reciprocal  Rela- 
tions. 

Apart  from  the  proposition  of  Onsager's  Reciprocal  Rela- 
tions, the  thermodynamics  of  irreversible  processes  aims  at 
establishing  the  true  "forces"  behind  the  processes.  The  correct 
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choice  of  conjugated  forces  and  fluxes  Is  itself  a  major  task  in 
the  application  of  irreversible  thermodynamics. 

2.  Diffusional  Phenomena:  In  diffusional  phenomena,  the 
classical  Pick's  law  attributes  the  cause  of  diffusion  to  the 
concentration  gradient.  The  diffusional  coefficients  thus 
defined  are  found  to  depend  on  the  concentration  and,  hence,  are 
not  true  constants.   In  the  formulation  of  irreversible  thermo- 
dynamics the  flux  is  written  as  linear  combination  of  certain 
forces  (as  in  Equation  (2-1)): 


n 

z 

3-1 


"i  =  .?,  LijXj  (2-1) 


inhere   the   Lj,»s   are  true   coefficients    in  the  sense   that  they  are 
independent   of  the  forces.     For   isothermal  diffusion  it  has  been 
proposed  that  the   gradient  of  chemical  potential  is   the   true 
cause   of  diffusion   (12),    (13).     Thus,   for  an  n-component  system 
under  isothermal  condition,    one  may  write: 

n 
Ji  =    •$!  LUV/V  i  =  1,    2,    .    .    .,    n  (3-1) 

where v/*     denotes  the  chemical  potential  of  component  j. 

The  present  chapter  is  concerned  with  the  question  of  the 
real  force  in  the  diffusional  processes.  If  the  diffusional 
coefficients  defined  in  terms  of  Equation  (1)  prove  to  be  true 
constants,  the  gradient  of  chemical  potential  should  be  con- 
sidered as  the  real  cause  of  diffusional  phenomena,  and  Equation 
(3-D  should  be  preferred  to  Fick»s  law.  Besides,  the  Onsager 
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Reciprocal  Relations: 


LLj  =  LjL  (2-5) 

could  be  used  with  benefit  to  reduce  the  experimental  work  In 
determining  diffuslonal  coefficients. 

Isothermal  diffuslonal  processes  are  described  mathemati- 
cally by  flow  equations  in  which  the  flux  or  current,  J.  of  the 
component  I  is  expressed  as  a  linear  combination  of  pertinent 
forces.   In  experimental  work  the  flux  J.  of  the  component  I  is 
equated  to  the  sum  of  concentration  gradients  each  multiplied  by 
a  diffusion  coefficient  DL .  of  one  type.  In  the  formalism  of 
irreversible  thermodynamics  the  flux  J.  of  the  component  I  is 
expressed  as  the  sum  of  the  gradient  of  chemical  potential  of 
each  component  multiplied  by  a  diffusion  coefficient  . .  of 
another  type.   In  the  following  the  Dj,  will  be  referred  to  as 
practical  diffusion  coefficients  and  the  . .  as  fundamental 
diffusion  coefficients.  Since  the  flux  J.,  may  be  defined  In 
various  ways  by  choosing  different  frames  of  reference,  the 
values  of  both  practical  and  fundamental  coefficients  depend  on 
the  frame  of  reference  selected.  Following  Fujita  (15)  and 
KIrkwood,  et  al.  (16),  one  Identifies  four  different  types  of 
frames  of  reference  capable  of  clear  definition.   In  the  defini- 
tions of  the  various  reference  frames,  one  dimensional  diffusion 
will  be  considered. 

3.  Frames  of  Reference:   In  defining  the  cell-fixed  frame, 
one  visualizes  a  cross-sectional  plane  1  cm2  In  area  and  per- 
pendicular to  and  fixed  with  respect  to  the  wall  of  the  cell. 
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The  flux   (Jj_)c    of  component   i   Ls   defined  as   the  number   of  grams 
of   L   crossing  the  mentioned  plane    in  one  second.      If  the  mean 
linear  velocity   of  the  component   i  x^ith  respect  to  the   cell   is 
denoted  by    (Uj_)c  and   the  concentration  in  gram  per  cc  by  C,,    it 
follows   from  the   definition  that: 

<°Vo  =   <VcCi  (3-2) 

One  may  choose  as  the  reference  plane  a  plane  moving  with  the 

mean  velocity  (Uq)c  of  the  solvent.  The  subscript  0  will  be 

used  to  represent  the  solvent  hereafter.  The  flux  (J.)0  of 

component  i  on  this  solvent-fixed  frame  is  related  to  (J.)  as 

i'  c 

follows : 

(Ji>0  =  (Jl>c  "  (VcCi  (3-3) 

By  substituting  i  =  0  into  Equation  (3-3)  and  noting  that 
(JQ)0  =  0,  one  obtains: 

tU0>c  =  (Vo/00'  C3-W 

and  Equation   (3)    takes   the   form: 

(jl)0  =  (jt)0  -  «yc0)(j0)c  (3.5) 

The  velocity  of  the  center  of  mass  relative  to  the  cell  is 

defined  as : 


W»  =  JbWo/g  °1  •  (3-6) 
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Making  use    of  Equation   (3-2)    and   the  expression  for  the   density 

n-1 
of  the   system  P=  X     Cj_,    one   rewrites  Equation   (3-6)    in  the 

following  form: 


(U)K0  =  <7J   1L    (Vc  07) 


One  now  defines  a  plane  moving  with  the  velocity  (U)n  as  the 
mass -centered  frame.  The  flux  (J^)™,  of  component  i  based  on 
the  mass-centered  frame  is  then  written  as: 

(JL)M=  (JjL)c  -  (U)MCCL  (3-8) 

Or,  substituting  Equation  (3-7)  into  Equation  (3-8),  one  has 

(Ji'H  -  (Vo  "  ( V*>  %   (Vc  (3-9) 

Finally  one  may  choose  as  the  reference  frame  a  plane  moving 
with  the  volume-centered  velocity  (IT)   which  is  defined  as: 

VC 

(U)vc  =  |0  7^liVOo/%  \°i  (3-10) 

where  7^  cc/s  denotes  the  partial  specific  volume  of  component  1. 
Since: 

n-1 

i?o  ^  =  l'  (3-n^ 

Equation  (3-11)  can  be  rewritten  as: 
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(u)vo  ■  lo  Wc 

(3-12) 

and 

the  flux,    (J»)v  °f  component   i  base 

d    on  the  volume 

fixed 

frame    is   in  the  form: 

(Vv  =    (V0    -    ^vo0!  ■    <Vc   " 

n-1 

c     X   V.(j  ) 

1   j=o      J      J   c 

(3-13) 

I}..     Relationships  Among  Fluxes  Bas 

ed   on  Various  Frames    of 

Reference:     The  relationships   among  the 

(JjOq,    UjOji*   and    (J  ) 

may 

be  derived  by  simple  manipulations. 

For  example,   combining 

Equations    (3-5),    (3-11),   and   (3-13)    one 

obtains : 

<jl'v  =  <Vo  -  ci  %Wo   . 

(Jo}o  =  ° 

(3-11+) 

i  =   1     2 

n-1 

The 

flow  equations  based    on  the   various 

frames   of  reference  are 

now 

written  as   follows: 

n-1 

.    .    . ,   n-1 

(1-15) 

n-1 

\ 

(1-16) 

-(Jt)v  =  J.'    (Dlj)y  C. 

(1-17) 

Here 

the  n(J\)    are   linearly  dependent,   hence   only   (n-1) 

fluxes 

appear  in  the   equation  and   the   negative 

sign  Is  adopted 

In 
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conformity  with  the   convention  that  the  positive  direction  of 
flux  is   along  the  direction  of  decreasing  concentration.     When 
Equations    (3-15)    and    (3-17)    are  substituted   into  Equation   (3-D|), 
one   obtains : 

<Vv  =   (D1350  "  Ci  £  Wo  »-18) 

In  Equation  (3-18),    the  diffusion  coefficients    (D.,  )_.  are 
expressed   in  terms   of  the   coefficients    (Dtj_)0.      One  now  seeks 
to  invert  Equation   (3-18)    in   order  to  have    (D, .)n  expressed   in 
terms    of   (*>..,)_.      To  this   end,    one  multiplies  Equation   (18)   by 
V.    and  sums   over  the  variable   i: 

&  71<Vy  =  "^x  VD1J>0  -  li  °1VL  4  Vk<Dtk>0  (3-W) 

Since    one  has   from  Equation   (3-5) : 

Equation  (3-19)   becomes: 


or 


Ix  Tl(D^0 :• o^  4  VVT  (3-20, 
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ReplacjLnc  the  dummy  Index  k  on  the  right  hand  side  of  Equations 
(18)  and  (20)  by  i  and  comparing  the  resulting  equations  one 
obtains  the  desired  relationship: 

(V<>  "  (DLj'v  +  0^  1J  7l(Dlj'v  (3-21) 

• 

5.   Legitimate  Choice  of  Conjugated  Forces  and  Fluxes:  At 
this  point  the  question  arises  as  to  which  set  of  fluxes,  when 
conjugated  with  the  gradient  of  chemical  potential,  will  give 
the  desired  or  the  legitimate  choice.  Since  the  invariance  of 
entropy  production  cannot  be  relied  on  In  the  determination  of 
the  correct  set,  either  experimental  verification  or  theoretical 
Investigation  of  microscopic  nature  will  serve  as  the  basis. 
Hooyraan  (17)  has  shown  that  at  thermal  and  mechanical  equili- 
brium the  solvent-fixed  frame  (Jj)0  is  the  desired  set  of 
f luxes . 

ITow,  among  the  various  set  of  fluxes,  the  cell-fixed  frame 
is  the  only  one  that  will  conveniently  yield  to  practical 
measurement.  However,  it  can  be  shown  (15)  than,  when  average 
values  of  partial  specific  volumes  7j_  can  be  used  with 
reasonable  accuracy,  the  volume -centered  frame  coincides  with 
the  cell-fixed  frame. 

When  one  writes  the  flow  equation  in  its  fundamental  form 
based  on  the  solvent-fixed  frame  of  reference: 

n-1 
:fJl}0  =  &  (iVoV/V     I  =  1,  2,  .  .  .,  n-l      (3-22) 
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The  resulting  (H.  )   are  the  coefficients  expected  to  obey 
Onsagerts  Reciprocal  Relations.   Moreover,  Equation  (3-22)  would 
be  considered  as  the  fundamental  flow  equation  giving  account  of 
the  real  "driving  force". 

6.  Gradients  of  Chemical  Potential  as  Driving  Forces:  It 
is  the  specific  purpose  of  this  chapter  to  examine  the  coeffi- 
cients (O     )      in  Equation  (3-22).  Previous  investigations  (13) 
have  been  confined  to  the  validity  of  Onsager»s  relations  among 
the  (A  .)  for  L  j£   j.  No  attention  has  thus  far  been  paid  to 
the  dependency  of  Cflj<)0  on  the  chemical  potential.  From  the 

practical  point  of  view,  the  dependency  of  (CI.  .)   on  the  cherai- 

*•  J    0 

cal  potential   is   a  more    important  problem  than  Onsager's 

Reciprocal  Relations  governing  the   cross   coefficients    (CI.  .) 

LJ'O 

for  L  f£   J,  which  are  of  less  significance  than  the  main 
coefficients  C^^q.  If  the  (n...)0  are  found  to  be  independent 
of  the  chemical  potential,  Equation  (22)  should  replace  Fick's 
law  as  the  fundamental  equation  of  diffusion. 

Expressions  for  the'  (^Lj-)0  in  terms  of  the  measurable 
(DLj)0  are  derived  as  follows.   One  first  rewrites  Equation  (22) 
by  expressing  ".  in  terms  of  the  gradients  of  concentration  as: 

Replacing  the  dummy  index  j  in  Equation  (3-15)  by  k  and 
comparing  the  resulting  .  express  ion  with  Equation  (3-23)  one 
obtains : 
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^Wo  ■   ?,    (i}ij}0a4f            i  =  1,    2,    .    .    .,    ft-l            (3-21,) 
J— J-                     k 

The    Inverted   equations   are: 

n-1 

(GLJ0  =   Z     BJk   (Dlk)Q,             L,    J   =   1,    2,    .    .    .,    n-1      (3-2^) 

where  E.n    are  elements    of  the    Inverse  matrixfBjof  the  matrix 

[  £l!  ]  Ln  which  JLl  is  an  element   of  the    1th  row  and  K  n  column: 
1  2  c  J                      ^  Ck 

iBj-t^r1 

and 

Bjk       * 

k,1 

(3-26) 

In  Equation   (3-26), 

dA 
dc 

is   the  determinant   of  the  matrix  1^.] 

Ld  c  * 

and 

dc 

is   the   cof actor  of  the  element  t— . 
kj                                                                   ^ 

For  three   component  system,   Equation   (3-25)    takes   the  forms: 

'      «%»o  ■  1 

5^11)0-^^12)0 

» 

(%2>0  =  Z 

• 

-zq  (Di2)0  -.rt|  (Dn}o 

$ 

<%>„  =  J 

'  ¥2           ¥2 

3c£  (D2l)o  "  JcT  (D22)q 

(3-27) 

**                                                                                                               * 

/ 

- 

% 
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<n22>0  =  X 


(Dp9)n   " 


ty 


22'0  ~a  c. 


(Dpi) 


21'0    » 


where 


A  = 


<W    «¥2        a/V  ^2 


The    Isothermal  diffusion  coefficients   have  been  experi- 
mentally determined  for  some   ternary  systems.     However,   there   is 
only  one  system,    Ha  CI  -  KC1  -  water  at  25°C,   for  which  complete 
thermodynamic  data  are  available  for  direct  computation  of  the 
tflLj)0  from  the    (Dj,    )Q. 

The  data  used   in  the   following  calculations   are   those 
summarized  by  Pitts    (9). 

For  ternary  systems,  Equations    (3-21)   are   in  the  form: 

c, 

V1(DU)V  +  V2(D21)y 


1 


(Dll>o  =   (Dll)v  +  ^T 

0    o 

C 


<D12>o  "   CD12)T  +  -Jr^CD^  +  V2(D22)    1 

•] 


(D21)0  =    (D01)__  + 


21/v  T  C0Y0 


(3-28) 


7l<Dll>v  +  VD21> 


(D00)      =    (D00)„  + 


22;o 


22yv 


o  ot  J 


Equations  (3-27)  and  (3-28)  are  used  to  compute  the  (fl.  .) 

iy  o 

from  the  data  of  Table  3-1.      IBM-1620  was  used  for  this,   and 
the  calculated  results   are  presented    in  Table   3-2. 

Prom  the  calculated  results   the  following  conclusions  may 
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Table  3-1.     Data  for  the  system  Na  CI  - 

KCl  -  HpO 

at  25°C 

for  four  different  compositions. 

* 

A 

B 

C 

D 

C0(g/cc)                                     0.9855 

0.9782 

0.9808 

0.9731+ 

C1(g/cc)                                     0.01I[61 

0.011+61 

0.02922 

0.02922 

C2(g/cc)                                    0.01863 

0.03727 

0.01863 

0.03727 

^0(cc/g)                                    1.002 

1.002 

1.0020 

1.002 

^(cc/g)                                    0.3197 

0.3270 

0.3273 

0.3339 

V2(cc/g)                                    0.3672 

0.3933 

0.391+1 

0.1:010 

(^-i)/RT                                 1.686 

1.536 

0.9502 

0.871+9 

{J-±)/m                          0.3791 

a  Cp 

0.2637 

0.2628 

0.20li.6 

(J£i)/*T                                 0.371(5 

0.2590 

0.2579 

0.1997 

<5^-)At                           1.007 

0.5586 

0.9176 

0.5129 

(Dll'v  *  lo5  (#5§>            l-37 

1.35 

1.1+1 

1.38 

(Di2»v  x  lo5  (m>       -°-00 

-0.00 

-0.00 

-0.00 

(D2l'v  x  lo5  <t^>            0.18 

0.21+ 

0.11 

0.20 

(D22>v  *  1qS    (^               1-82 

1.83 

1.81 

1.83 

-"-The  subscripts   0,    1,   and   2 

refer  to  Hp0 

,   Na  CI,   and  KCl 

respectively. 

% 
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Tab 


le  3-2.     Values    of   (fl.  .)      for  the   system  Fa  CI  -  KC1  - 


HO  at  25°C   for  four  different  compositions. 


A 

B 

c 

D 

CQ(s/cc) 

0.98^5 

0.9782 

0.9808 

0.973)| 

Cj^Cg/cc) 

O.OIJ46I 

0.011*61 

0.02922 

0.02922 

C2(g/cc) 

1 

0.01863 

0.03727 

0.01863 

0.03727 

t°U>o* 

10* 

•   RT 

0.8890 

0.9570 

1.619 

1.71:4 

(n12)0:c 

10S 

*   RT 

-0.32k3 

-0.1i325 

-O.J4I1.O6 

-0.6531 

(n21)o  x 

10* 

r  RT 

-0.3185 

-0.1}  21^9 

-O.J4lj.83 

-0.5217 

(O>2)0  x 

IflS 

•   RT 

1.9i|0 

3.525 

2.115 

U.323 

28 
be  drawn.  While  the  Onsager  Reciprocal  Relations:   Ul.p)  = 
(.0  . )  ,  are  verifiable,  the  values  of  the  (fl.   . )   are  found  to 
vary  with  the  composition  of  the  system*  The  significant  fact 
is  that  the  dependency  of  the  (Q.  .)      on   the  composition  is  even 
more  marked  than  that  of  the  conventional  Pick's  law  constants 
(Dj,)  .  The  unexpected  results  may  be  ascribed  in  part  to  the 
inaccuracy  in  the  chemical  potential  data  involved  in  the  calcu- 
lation. However,  the  variance  of  the  so-called  fundamental 
coefficients  is  so  significant  that  there  seems  to  be  ample  room 
for  doubt  as  to  the  feasibility  of  replacing  Pick's  law  by 
Equation  (3-22)  as  the  fundamental  equation  for  diffusion. 

7.  Summary:  The  attempt  to  establish  a  fundamental  equa- 
tion for  isothermal  diffusion  based  on  the  scheme  of  irreversi- 
ble thermodynamics  has  yielded  negative  results.  As  a  by- 
product of  this  'attempt,  however,  Onsager's  Reciprocal  Relations : 

^ij^o  =  ^ji^o*  Ln  thLs  case>  arG  found  to  hold  reasonably  well 
in  agreement  with  Miller's  work  (13). 
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Chapter  1;. 

A  FORMALISM  DUE  TO  TYKODI  AMD   ITS  APPLICATIONS 

1.      Introduction:      In  the  present  chapter  a  new  approach 
developed   by  R.    J.   TykodL   is   presented.      The   original   treatises 
have  appeared   in  a  series    of  papers    (18),    (19)    in  the  Journal 
of  Chemical  Physics.      In  dealing  with  non-equilibrium  phenomena 
TykodL1 s   method  aims   at  preserving  the  phenomenological  spirit 
characteristic   of  the  classical  equilibrium  thermodynamics.     The 
physical  quantities    involved   in  this  formalism  are    introduced   in 
as    operational  a  manner  as  possible.     The   scheme   of  presentation 
in  this   chapter  will  follow  TykodL* s  manuscript   (20)    of  a  forth- 
coming book.      It   is   essential  that  the   closed  system  approach  is 
adopted  to  formulate  a  problem  by  TykodL* s  method.      In  contrast 
to  the   open-system  method  which  will  usually  provide  a.  more  clear 
cut  approach  to  typical  engineering  problems,   the   closed-system 
method  may  tax  the   imagination  more   in  defining  a  system.      In 
formulating  a  problem  the  system  is  defined   in  such  a  way  as   to 
ensure  that   only  heat  and  work  can  be  exchanged  between  the 
system  and   the  surround ines . 

2.     Basic  Technique:     The  basic   technique  employed   in 
analyzing  a  problem  by  TykodL* s  approach  Ls   to  define  a  system 
in  such  a  way  that  the  system  may  be   separated   into  terminal 
parts   and  parts    on-the-line    (betx^een  terminal  parts).     The 
thermodynamic  properties   are  uniform  throughout  each  terminal 
part,   and  the  gradients    of  the   thermodynamic  properties   exist 
only   in  the  parts    on-the-line.     The  parts    on-the-line  are   taken 
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to  be  point-wise  time  Invariant.  In  this  respect  a  severe 
restriction  is  imposed  upon  Tykodi's  method.  This  approach  is 
confined  to  the  treatment  of  steady  rate  processes.   In  each 
case  a  definite  amount  of  mass  is  taken  as  the  system  (See 
Pig.  1,  with  accompanying  illustrations). 

Pour  basic  assumptions  are  made.  They  are: 

1.  The  system  is  capable  of  sustaining  steady-rate  opera- 
tions under  the  described  conditions. 

2.  The  system  can  always  be  divided  into  terminal  parts 
and  parts  on-the-line,  the  parts  on-the-line  being  stationary 
during  steady-flow  operations. 

3.  Certain  non-flow  states  can  be  treated  as  being  the 
limit  of  sequence  of  steady-flow  operations. 

l[.  In  a  sequence  of  steady-flow  states  converging  on  a 
non-flow  state,  the  limiting  state  is  a  state  of  minimum  entropy 
production  relative  to  the  neighboring  flow  states. 

One  may  interpret  the  preceeding  set  of  assumptions  in  the 
following  manner.  The  first  three  assumptions  are  pure  assump- 
tions made  on  the  physical  situation.  These  three  assumptions 
act  together  to  impose  an  inherent  limitation  upon  the  scope  of 
applicability  of  this  method.   The  fourth  assumption  is  of  a 
rather  different  nature.   It  should  be  considered  as  a  fundamen- 
tal postulate  upon  which  the  validity  of  the  whole  scheme  of 
Tykodl»s  approach  is  hinged.  Tykodi  has  suggested  (20)  that 
the  proof  of  assumption  Ij.  may  be  obtained  through  statistical 
mechanical  investigations.  As  a  matter  of  fact,  assumption  i| 
is  known  to  hold  in  monothermal  fields  (i.e.,  the  system  is  in 
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Fig.  I.  A  mass  of  gas  in  adiabatia  flow 
across  a  porous  plug,  P+^  P-  . 
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contact  with  a  single  heat  reservoir)  when  the  reference  non-flow 
state  is  in  a  state  of  thermodynamic  equilibrium;  since  the  rate 
of  entropy  production  is  zero  in  a  state  of  thermodynamic 
equilibrium,  whereas,  in  any  neighboring  steady-flow  state,  the 
rate  of  entropy  production  must  be  positive  definite. 

Tykodi's  theory  of  non-equilibrium  phenomena  may  be  said  to 
be  an  exploration  of  the  consequences  of  assumption  h   for  gen- 
eral cases.   In  this  chapter  Tykodi»s  fourth  assumption  will  be 
referred  to  as  the  fundamental  postulate. 

3.  Thermodynamic  Relations:  In  the  following  treatment 
the  standard  notations  for  thermodynamic  functions  in  classical 
equilibrium  thermodynamics  are  retained.  Thus,  U,  H,  A,  G-,  and 
S  will  denote  respectively  the  total  energy,  the  enthalpy,  the 
Helmholtz  free  energy,  the  Gibbs  free  energy,  and  the  entropy. 
Q  and  W  will  represent  heat  and  work  supplied  to  the  system.  A 
superior  dot  will  be  used  to  denote  the  time  derivative  or  rate. 

XT'  (system)  =  Q  +  W  (lj.-l) 

A  distinction  is  drawn  between  monothermal  and  polythermal 
processes  (and  systems).  A  monothermal  process  is  a  process  in 
which  the  entire  system  is  in  heat  communication  with  a  single 
heat  reservoir  of  temperature  T  during  the  process.  A  polyther- 
mal process  is  a  process  in  which  the  terminal  parts  of  the 
system  are  separately  in  heat  communication  with  heat  reservoirs 
of  temperature  Tj.  during  the  process.   Heat  exchange  between  the 
surroundings  and  the  parts  on-the-line  may  be  permissible  but  is 
not  analyzed  further  in  the  present  work.   For  a  monothermal 
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process,  then: 

4  -  -T  S(r)  (h-2) 

(v) 
where  Sv  '   denotes  the  rate  of  accumulation  of  entropy  in  the 

heat  reservoir.  Similarly,  for  a  polythermal  process 

\   =  -Zt.S.W  ([,_3) 


Q 


i  &  *■ 


When  a  steady  exchange  of  mass  between  terminal  parts  of  the 
system  takes  place,  the  parts  on-the-line  is  point-wise  invariant 
in  time  during  the  steady-flow  process.  Under  such  conditions 
one  writes : 

U  (system)  =  Z  Ufa  (l±-l±) 

1 

where  U  represents  the  total  energy  per  mole  of  the  material  in 
the  i   terminal  parts  of  the  system. 

As  an  illustrative  example  the  classical  porous  plug  experi- 
ment is  analyzed  here  by  the  present  method  (21).  Pig.  1  depicts 
the  classical  Joule-Thomson  experiment.  A  definite  amount  of 
gas  and  a  plug  constitute  the  system  and  are  thermally  insulated 
from  the  surroundings.  The  usual  assumptions  are  made  that  the 
gas  is  flowing  across  the  porous  plug  at  a  steady  rate  of  N  moles 
per  second,  that  the  gradients  in  temperature  and  pressure  can 
be  localized  in  a  region  in  the  close  vicinity  of  the  porous 
plug,  and,  farm  from  the  plug,  the  thermodynamic  state  of  the 
gas  is  characterized  by  the  variables  (T+,  P+)  and  T_,  P_). 
Under  such  conditions  one  writes: 


• 

U   (system) 

=  W  =   -P . V.    +   P  V 

+    +             -    - 

3k 

(li--5) 

Since  the   flow 

is  between 

terminal 

parts   and 

the 

part   on- 

-the-line 

is 

point-wise   t 

* 

U   (system) 

ime    invariant,    one 

=  ttJt  +  tj  n 

+    +              —    — 

has 

! : 

(lj.-6) 

where   the  superior  bar  inc 

iicates   a 

molar  property 

'  of  the 

gas. 

the 

■■  steady  flow 

-N     =  Jl    = 

+ 

condition 

11 

dictates 

• 
• 

(k-7) 

Hence   one   obtains 

ro_  -  rr+)  i 

I  =  (P+V+  - 

■  PJM   N 

Ui--6) 

or 

ff_  =  ff+ 

■  &-9) 

whe 

re 

H  =  U  +  PV 

• 

In 

line  with  Equation   (1{-L\ 

.)    one  can, 

in  general, 

write 

• 

U  (system) 

-TS^   + 

W 

Ut-10) 

• 

H   (system) 

=  |hLlTL   = 

-TS^   + 

w 

+      Tiii 

Ol-ll) 

• 

A    (system) 

=fVTi  ■ 

-T  (£W 

+ 

S^)    + 

• 

W 

(U-12) 

• 

G    (system) 

-T    (S(^ 

+ 

sLNL)  + 

w  + 

PA 

(U-13) 

3S 

The  rate  of  entropy  production  is  given  under  steady  flow 

conditions  by: 

• 

G  =  S(r)  +  SjNj,          \       •                    (k-Th) 

• 

In  the  following  treatments  It  will  be  assumed  that,  for  all  the 

•            • 

systems  to  be  considered,  W  +  Pj.Vj_=  0,  i.e.,  the  only  work 

involved  is  the  PV  work.  Under  such  circumstances  Equation 

(It- 13)  takes  the  form: 

-G  (system)  =  TCT                                  (U-15) 

It.   Migratlonal  Equilibrium:   The  conditions  for  a  mlgra- 

tlonal  equilibrium  will  now  be  sought.  By  migratlonal  equili- 

brium one  understands  a  steady  state  corresponding  to  a  non-flow 

state  with  respect  to  a  given  chemical  substance.   When  a  steady 

migration  of  a  given  chemical  substance  Into  the  K   terminal 

part  of  the  system  takes  place,  one  writes: 

-G  (Hk)  =  Ttf(Nk)                                  (Ii-16) 

For  the  non-flow  state  with  respect  to  the  given  chemical  sub- 

stance It  follows  that: 

-G(0)  =  Ttf(0)                                      (14.-17) 

Combining  Equations  (16)  and  (17)  one  has: 

-Ag  =  T^<r                            (i[-i8) 

- 

« 
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where 

AZ  =  Z  (\)    -  Z  (0) 

Dividing  both  sides  of  Equation  (I4.-I8)  by  N,  and  passing  to  the 
limit  of  vanishing  mass  flow  via  a  sequence  of  steady-flow 
states  one  writes: 

« 

-    .lira     (AG/4)   =  T    .lira     (A<V?U  (^-19) 

K-"  0  K  K— 0 

k  k 

Since   this   notion  of  limiting   operation  is  a  basic   tool  in  the 
application  of  Tykodi's-  formalism,    one  wishes   to  employ  a  short- 
hand  notation.      One  adopts   the  following  notation: 

6Z/S\  =    .Um    j  (z(Nk)    -  Z(0)]    Ak|  (l|-20) 

k        L 

the  limit  being  always  approached  via  a  sequence  of  steady-flow 
states.  The  fundamental  postulate  now  takes  the  form: 

**/l\   =  0  .  (11.-21) 

Substituting  Equation   (1+-21)    into  Equation   (I4.-I9)    one    obtains: 

-&i/6\  =  -SgtNt)./Sk  =  0  (Li-22) 

For  the  case  of  a  steady  flux  of  mass  between  two  terminal  parts 
{ot,/S)    of  a  system,  one  has  N*  +  IT^  =  0  and  that  Equation  (4-22) 
reduces  to: 

8*-*,-°  U1--23) 

Equation  (23)  yields  familiar  results  in  equilibrium  thermodyna- 
mics. For  example,  if   and   refer  to  gas  phase  and  liquid 
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phase   respectively,    the  condition  for  the  monothermal  gas-liquid 
equilibrium  is    just: 

gas  liq 

■where/*  stands  for  the  chemical  potential,  or  the  partial  molar 
Gibb's  free  energy. 

5.  Migrational  Equilibrium  in  a  Monothermal  Concentration 
Field:  The  problem  of  migrational  equilibrium  in  a  monothermal 
concentration  field  will  now  be  analyzed.  The  situation  is 
depicted  in  Pig.  2. 

There  are  two  fluid  components  1  and  2  present  in  the 
spatial  region  marked  MIX,  which  communicates  with  terminal  parts 
A  and  E  via  membranes  permeable  to  component  1  alone  and  with 
terminal  parts  oL  and  j3   via  membranes  permeable  to  component  2 
exclusively.  A  steady  flow  of  component  1  is  induced  through 
the  MIX  region,  and  the  conditions  for  the  migrational  equili- 
brium of  component  2  in  the  concentration  field  formed  by  the 
flow  of  component  1  are  sought. 

Now  with  respect  to  the  fundamental  postulate,  the  refer- 
ence state  of  the  system  corresponds  to  the  nonflow  state  of 
component  2  (with  the  flow  of  component  1  alone,  or  monoflow 
state).  When  a  steady  flux  of  component  2  is  established 
through  the  MIX  region  with  the  pressures  PA  and  P  kept  fixed, 
one  writes  for  the  biflow  state: 

-C-(%,N<)  =  T  (1^,^)  (1^-25) 
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Fig. 2.  Migrational  equilibrium  in  a  mono- 
thermal  concentration  field. 
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At  static  condition,  the  partial  molar  G-ibbs  free  energy  of  a 
chemical  substance  is  identical  to  the  chemical  potential.  When 
a  substance  is  in  motion,  one  adds  the  molar  kinetic  energy  of 
macroscopic  motion  K,  to  the  chemical  potential  to  get  the 
partial  molar  C-ibbs  free  energy  (21)  and  one  writes: 

g  =/*+  K  (U-26) 

An  important  property  of  the  molar  kinetic  energy  of 

macroscopic  motion  is  obtained  as  follows.  Allowing  a  bracketed 

subscript  (i)  to  denote  the  property  of  component  i  and  denoting 

the  molecular  weight  and  volumetric  flow  rate  by  M(j\  and  V/;\ 

respectively  and  if  the  cross  sectional  area  of  the  conduit  . 

perpendicular  to  the  flow  direction  is  denoted  by  A   one  has: 

3     - 

M  v" 

(1)    /(1).2 

Hi)  =-?-  ncro 

s 
or 


K(i)   =^-  (N(i)V(1)/As)2 

•        2         -  &-271. 

2  "fW   a      ; 

s 


and 


huKu)  =  %)3  <f>  4^> 


H 


Ij-o 

Applying  the   limiting  flow  operation: 

^(i)K(i)Mk)    =32T(j.)2    (§)     (-P-)2    (^(l)/S(k)) 

s 

=  3K(1)  (^n(1)/Jn(1c)5  U:--28) 

The   left  hand  side   of  Equation   (Li-25)    is,    in  a  more   explicit 

form, 

Because  of  the  nature  of  the  problem,  the  E  terms  for  com- 
ponent 2  may  be  neglected.  And  one  rewrites  Equation  (2j.-2f>)  as: 

•  •  •  • 

How,   noting  that  NA  +  Kg  =  0   and  IT    +  N     =0  and   applying  the 

limiting   operation,   °/6  HUf    to  Equation   (lf-29),    and  using 
Equation   (ij-28)    one  has: 

-  ( H  +  3V  ■  ^b  +  3V)  <ty'V  +  <A  ■  >*>}  = T  £ 

(1l-30) 
The  right  hand  side  of  Equation  (30)  vanishes  according  to  the 
fundamental  postulate.  Hence,  one  has: 


hi 

where 

g*  =  /*  +  3K  and  (  R^J  =  £%/?£, 

Identifying  subscripts  are  attached  to  the  (  rJ  quantity  to 
Indicate  the  points  in  the  concentration  field  between  which 
chemical  potential  differences  are  being  computed.  If  the  flow 
of  the  two  fluids  is  coupled,  i.e.,  (R^aJ  ^   0,  then  the  flow  of 
component  1  sets  up  a  gradient  in  the  chemical  potential  of  com- 
ponent 2.  Here  g  behaves  as  a  driving  potential  for  the 
exchange  of  mass.  The  quantity  g  =M  +  3K  plays  an  important 
role  in  subsequent  applications. 

6.  A  Generalized  Problem  of  Migrational  Equilibrium:  A 
situation  as  illustrated  in  Pig.  3  will  be  analyzed.  Pig.  3 
depicts  a  one-component  fluid  system  consisting  of  two  parallel 
streams  of  fluid  (floxtfing  through  conduits  of  uniform  bore) 
connected  by  a  linkage  through  which  a  steady  flux  of  fluid 
takes  place.  The  system  is  divided  up  into  terminal  parts 
(x+,  x-)  with  well  defined  fluid  velocities  under  all  conditions 
and  parts  on-the-line  which  are  stationary  under  conditions  of 
steady  flow,  including  conditions  of  steady  exchange  of  fluid 
between  the  tubes  through  the  connecting' linkage.  Por  this 
situation,  Equation  (li-lk)  is  substituted  into  Equation  (I|-13) 

-Is^  =  T<T  -  W-lP^    x  «  *,  f,,  7  B  +,  -    ([^33) 

where  g       =  U     +  k       and   $  represents   the   entropy  production  of 
the  process. 
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Fig.  3.  Analysis  of  a  general  problem  in 
migrational  equilibrium  .PX+^PX-, 
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In  the  presence  of  a  steady  flux  of  fluid  through  the 
connecting  link,  one  writes: 

**  +  K-  -  **  <    i +  % =  °  .  (WU) 

where  V  represents  the  rate  of  Influx  of  fluid  Into  the  x  con- 

dult  via  the  connecting  link.  Equation  (l|-33)  is  now  written  in 
an  expanded  form: 

Applying  the  limiting  operation,  6/£\t$    *©  Equation  (37),  and 
using  the  result  of  Equation  (28)  i.e.,  £(NX+Kx+  +  \J-xJ/^x   = 
3K  ,  (K„  being  considered  as  the  average  value  of  E^.  and  K^J 


x*  %  x 


where  g/  =/x  +  3KX  and  (dJ  -  -  8  C^VV^V  Equatlon  (36) 
is  the  final  result  of  the  preceeding  analysis.   It  has  a  far 
reaching  significance  in  the  analysis  of  migrational  equilibrium 
in  general.  The  [v]    term  may  be  named  as  "dissipation  function" 
and  physically  represents  the  dissipating  effect,  i.e.,  the 
viscous  dissipation. 

Now,  one  wishes  to  explore  the  implications  of  Equation  (36) 
on  the  problem  of  phase  equilibrium  in  the  most  general  situa- 
tion, i.e.,  the  various  phases  in  contact  are  allowed  to  move 
with  different  velocities.   One  considers  the  distributions  of 
a  certain  chemical  substance  between  two  phases  which  are  in 
relative  motion. 


Fig.  k   shows  the  migrational  equilibrium  (or  equilibrium 
distribution)  of  a  substance  between  two  phases  in  relative 
motions.  Since  the  interphase  of  the  two  phases,  as  pictured 
in  Pig.  k,  may  be  thought  of  as  an  "extended"  connecting  link 
(referring  to  Fig.  3),  one  expects  that  the  equilibrium  distri- 
bution will  be  governed  by  Equation  (36): 

tf-ftJ-i/-[i»>J-  '  (,l-36) 

If  one  is  allowed  to  neglect  the  dissipative  term  D  ,  one  has: 

for  the  condition  of  migrational  equilibrium. 

When  one  may  further  neglect  the  kinetic  effect  of  macro- 
scopic motion,  K,  the  familiar  relationship  established  by  the 
classical  equilibrium  thermodynamics  results: 

A=A*  (4-38) 

In  the  following  sections  use  will  be  made  of  the  proceeding 
results  in  analyzing  two  simplifying  limiting  cases. 

7.   Corrections  for  Manometer  Readings:  The  problem  of 
measuring  gas  flow  by  use  of  the  mercury  manometer  will  be 
analyzed.  Under  the  steady  flow  condition  the  gas  in  the  mano- 
meter arm  will  be  in  migrational  equilibrium  with  the  gas  flow- 
ing in  the  conduit.  Application  of  Equation  (37)  in  this 
situation  leads  to: 

/*  +  K  =/^man  U-39) 
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1+6 
inhere  unsubscripted  symbols  represent  the  mainstream  properties . 
When  the  state  of  the  gas  may  be  described  by  the  ideal  gas  law, 
one  has  (22) : 

JUL   =#°  +  RT  In  P  (1H+0) 

2 
and  K  =  i— — ,  where  H  and  v  represent  the  molecular  weight  and 

linear  velocity  of  the  gas,  Equation  (J|l)  yields: 

p/pman  =  exp  (-3Mv2/2RT)  (Il-4^1) 

To  examine  the  consequence  by  way  of  numerical  computation 
one  considers  the  case  of  air  flowing  with  a  velocity  of 
lOOm/sec  at  0°C.  Using  M  =  29  for  air: 

IIy2        3'x  29  x  (10, OOP)2 
2KT~  -  2  x  82.06  z  1033  k   9B0  x  273 

=  0.192if 
and  consequently 

P/P  „  =  e-°'19*     ■ 
'  man 

=  0.825 

In  the  proceeding  analysis  one  neglects  the  dissipative  effect 
(Dj,  and  thus  obtains  a  somewhat  exagerated  result.  However, 
this  example  may  serve  to  bring  out  the  effect  of  the  kinetic 
term  of  macroscopic  motion  when  the  gas  is  flowing  with  a 
velocity  of  the  order  of  lOOm/sec. 
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0.  Mass  Transfer  Processes:  In  the  usual  air  liquif ication 
process,  compressed  air  Ls  passed  through  packed  bed  dryers 
(usually  si.li.ca  or  aluminum  gels  and  more  recently  molecular 
sieves)  'prior  to  entering  the  refrigeration  unit.  The  velocity 
of  air  is  usually  high  enough  to  justify  consideration  of  the 
kinetic  term  K.  One  is  led  to  suspect  that  the  equilibrium 
moisture  content  of  the  flowing  air  would  be  lower  than  if  the 
air  is  in  static  contact  with  the  adsorbent.   One  pictures  the 
situation  by  means  of  a  model  in  which  the  adsorbent  is 
surrounded  by  a  thin  stagnant  film  of  air  beyond  which  the  bulk 
of  air  is  flox^ing.   Neglecting  the  dissipative  effect  one 
applies  Equation  (I)--37)  to  this  situation: 

where  the  unsubscripted  symbols  represent  the  mainstream  and  the 
subscript  s  refers  to  the  stagnant  air  film.  Under  such  condi- 
tions water  vapor  is  not  expected  to  obey  the  ideal  gas  behavior, 
and  one  xvrites  (22)  : 

/  ~/*°  +  RT  In  f  (lj.4,5) 

Tihere^tf      depends    on  the  temperature  alone,   and  f  denotes   the 
fugaclty. 

Substituting  of  Equation  (k~k$)    Into  Equation  (k-kk)   leads 
to 

RT   In  f     =  -  3K  =  -  m?  (U-k6) 


or 

f  =  exTD  I   -  3Mv2^ 
s 

f  is  the  saturation  fugacity  of  water  vapor  in  static  condition, 
s 

In  leaching  processes,  solvent  flows  past  stationary  beds 
of  solid  material.  Utilizing  the  similar  model  as  in  the  pro- 
ceeding example,  one  envisages  a  stagnant  solvent  film 
surrounding  the  solid  particle.  For  this  situation  one  writes 
(22): 

JU  aft*   +  RT  In  Y  x 

where//'"  is  a  function  of  temperature  and  pressure,  and  x  and 
stand  for  the  sole  fraction  and  activity  coefficients  of  the 
solute.  Carrying  through  similar  manipulations,  one  obtains: 

The  results  of  the  proceeding  analysis  point  to  a  reconsid- 
eration of  the  so-called  "driving  force"  in  mass  transfer.   In 
the  conventional  design  approach  to  the  design  of  mass  transfer 
processes  between  two  phases,  concentration  (or  partial  pressure) 
difference  in  the  two  phases  accounts  for  the  driving  force. 
The  equilibrium  concentrations  and  partial  pressures  are  tradi- 
tionally evaluated  at  static  conditions,  (i.e.,  when  there  is  no 
relative  motion  between  the  phases).  Equations  (U-I4.6)  and 
(k.-k7)   point  to  a  need  for  correction  on  the  values  of  equili- 
brium concentration  when  relative  motion  between  the  phases  is 
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significant.   It  Is  desirable  to  draw  a  distinction  between 

steady  state  concentration  (or  partial  pressure)  and  static 

equilibrium  concentration  (or  partial  pressure)  and  the  former 

should  be  used  whenever  relative  motion  between  different 

phases  is  of  considerable  significance. 

> 

9.  Summary:  Tykodi's  theory  has  been  successfully  applied 
in  the  analysis  of  migratlonal  equilibrium.  However,  it  must 
be  emphasized  that  the  legitimacy  of  Ty.kodi«s  theory  rests  upon 
a  postulate  which  is  yet  to  be  proven.   Moreover,  in  order  to 
apply  the  fundamental  postulate  in  the  form  of  the  limiting 
operator  *»—  it  is  essential  that  one  non-flow  state  with 

respect  to  a  certain  chemical  substance  must  be  identifiable. 
This  requirement  imposes  a  practical  limitation  upon  the  domain 
of  applicability  of  Tykodi's  formalism. 
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NOMENCLATURE 

A   =  Helmholtz  free  energy 

Ag  =  cross  sectional  area 

a   =  generalized  thermodynamic  function  (Chapter  1),  molar 

Helmholtz  free  energy  (Chapter  I|.)  ( 

B   =  magnetic  Induction 

c   =  concentration 

D   =  dissipation  function 

Dj_j  =  Pick's  law  dlffuslonal  coefficient 

f   =  fugaclty 

G   =  Glbbs  free  energy 

g   =  molar  Glbbs  free  energy 

H   =  enthalpy 

h   =  molar  enthalpy 

J.  =  flux  of  component  I 

K   =  kinetic  energy  of  macroscopic  motion 

Lt  .  =  constant 

M   =  molecular  weight 

• 

N   =  molar  rate  of  flow 

P  =  pressure 

Q  =  heat  supplied  to  a  system 

Pi  =  gas  constant 

S  -   entropy 

s  =  molar  entropy 

T  =  temperature 


U   =  wave  function  (Chapter  1),  mean  velocity  (Chapter  3), 

total  energy  (Chapter  I4.) 
u   =  molar  total  energy  \ 

V  =  volume  rate  of  flow 
IT   =  mean  molar  volume 

v   =  linear  velocity 

W   =  work  done  on  the  system 

x,   =  generalized  thermodynamic  force 

V  =  activity  coefficient 
/£       =   chemical  potential 

il.  .  =  fundamental  diffusional  coefficient 
(J)       -   angular  velocity 

Subscript 

C   =  cell  fixed  frame  of  reference 

0   =  solvent  fixed  frame  of  reference 

V  =  volume  centered  frame  of  reference 
M   =  mass  centered  frame  of  reference 
i,j  =  components 
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ABSTRACT 

This  work  consists  of  four  chapters.  In  the  first  chapter, 
the  context  and  scope  of  nonequi libr ium  thermodynamics  are  out- 
lined. Critical  considerations  are  given  to  the  criteria  of 
irreversibility.  A  brief  historical  development  of  the  subject 
is  sketched.   Onsager's  theory  and  associated  problems  are 
presented  from  a  practical  point  of  view.  The  second  chapter 
is  concerned  with  a  controversial  problem  on  the  application  of 
Onsager's  reciprocal  relations.   It  is  emphasized  that  a  legiti- 
mate choice  of  the  conjugated  forces  and  fluxes  depends  upon  a 
more  fundamental  criterion  than  the  notion  of  invariance  of 
entropy  production.  Chapter  three  deals  with  isothermal 
multicomponent  diffusion.  The  classical  Flck's  law  of  diffusion 
is  compared  with  a  formulation  based  on  the  scheme  of  nonequili- 
brium  thermodynamics.  The  system  NaCl  -  KC1  -  H20  is  investi- 
gated. While  numerical  computations  based  on  existing  data 
yield  favorable  results  for  Onsager'3  reciprocal  relations,  a 
vital  defect  in  the  formulation  based  on  the  scheme  of  non- 
equilibrium  thermodynamics  is  revealed.   In  the  final  chapter, 
a  nei*  approach  developed  by  R.  J.  Tykodi  is  presented.   In 
treating  nonequi libr ium  phenomena,  Tykodi' s  method  alms  at 
preserving  the  phenomenologlcal  spirit  characteristic  of 
classical  equilibrium  thermodynamics.  Problems  of  migrational 
equilibrium  are  successfully  attacked  by  Tykodi' s  method. 


